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ABSTRACT: 

An  analysis  for  the  transient  axial  response  of  a  gun  launched  motor 
case  with  unbonded  propellant  is  performed.  The  present  effort  extends  a 
previous  work  by  including  higher  harmonics  in  the  representation  of  the 
breech  force  at  the  aft  end  of  the  motor  case.  Including  the  higher  har¬ 
monics  leads  to  a  more  accurate  representation  of  the  actual  breech  force. 
The  results  of  the  analysis  indicate  that  the  dynamic  effects  on  the  in¬ 
ternal  axial  force  in  the  case  are  negligible  for  a  breech  force  duration 
of  0.026  seconds.  This  is  the  same  conclusion  that  was  made  in  the  pre¬ 
vious  work  based  upon  a  simple  representation  of  the  breech  force. 
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NOTATION 


area  of  motor  case  cross  section 
Fourier  cosine  coefficients  of  breech  force 
Fourier  sine  coefficients  of  breech  force 
speed  of  wave  propagation 
axial  strain 

Young's  modulus  of  elasticity 

axial  force 

constant 

length 

mass  of  motor  case 
breech  force 
time 

deformation 

displacement 

rigid  body  displacement 

transformation  variable 

transformation  variable 

coordinate  axes 

coordinate  axes 

time  duration  of  breech  force 

constant 

mass  density  of  the  motor  case  material 
axial  stress 


constant 


INTRODUCTION 


This  report  examines  the  transient  axial  response  of  a  gun-launched 
rocket  motor,  with  unbonded  propellant,  during  the  launch  phase.  It  is 
part  of  a  continuing  effort  to  develop  structural  design  procedures  for 
gun-launched  rockets.  In  a  previous  report  Uef.  1),  an  analysis  for  the 
time-dependent  axial  force  distribution  in  the  motor  case  was  performed 
in  which  the  force  at  the  breech  end  of  the  case  was  approximated  by  a 
one-half  sine  wave  over  the  duration  of  the  loading.  That  analysis  estab¬ 
lished  that  for  the  one-half  sine  wave  input,  a  time  duration  of  .012 
seconds  leads  to  negligible  dynamic  effects,  i.e.,  for  a  one  half  sine 
wave  of  .012  seconds  duration,  the  internal  axial  force  at  any  location 
along  the  case  appears  as  a  duplicate  of  the  breech  force,  but  of  smaller 
magnitude.  It  was  shown,  however,  that  decreasing  the  time  of  the  force 
duration  to  .0012  seconds  leads  to  significant  dynamic  effects.  This  re¬ 
sult  indicated  that  higher  harmonics  in  the  breech  force  could  cause  im¬ 
portant  dynamic  effects.  Since  the  actual  breech  force  is  lot  a  simple 
half  sine  wave,  but  in  fact  is  a  complicated  function  of  ti.ie,  an  investi¬ 
gation  was  made  to  consider  the  problem  further. 

To  account  for  the  higher  harmonics  in  the  breech-force,  the  analysis 
decomposes  the  actual  force  into  a  Fourier  series  with  27  terms.  This 
force  is  applied  to  the  case,  and  the  internal  axial  force  is  computed  at 
several  locations  along  the  case  over  the  duration  of  the  loading.  This 
analysis  then  establishes  the  contribution  of  the  higher  hermonics  in  the 
breech  force  to  the  dynamic  stresses  in  the  bar.  The  analysis  is  restricted 
to  the  free-free  motor  case. 
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ANALYSIS 


In  the  figure  below,  AB  is  the  motor  case,  treated  as  a  one¬ 
dimensional  bar,  under  the  applied  breech  force  P(t).  The  actual 
displacement  of  any  point  of  the  bar  with  respect  to  the  X,  Y, 
fixed  coordinate  system  is  denoted  ii.  This  motion  consists  of  two 
parts,  a  rigid  body  motion  u  and  a  strain-causing  deformation  u. 


♦  Y 


Then 


P(t) 


If  m  is  the  total  mass  of  the  bar,  then  the  rigid  body 
motion  is 

A  ,  EM  ... 

dt2  *  (2) 

The  differential  equation  of  motion  for  any  point  of  the  bar  is 

2  2 
3  G  r2  3  U  n 

^'c  ^-°  (3) 

where  c  is  the  speed  of  wave  propagation,  and  x  is  a  coordinate 
axis  parallel  to  x  moving  with  velocity  u.  The  constant  c  is 
given  by 

c  =  Vf  (4) 

where  E  is  Young's  modulus  of  elasticity  of  the  bar  material  and 
p  its  specific  gravity.  Combining  equations  (1)  and  (3)  gives 


(u  +  u)  -  c2  2-4  = 

at  •'-* 


ax 


2 


a2u 


Substituting  for  ,2  from  (2)  yields. 


at 

£ 

3t‘ 


u.  .  P(t)  .2  S  u 

O  '  _  “  ^  A  “ 


0. 


ax 


With 


2  2 

the  notation  (  )xx  denoting  4IJ.  ,  and  (  )  denoting  ^-4-4  , 

ax  at 


the  governing  differential  equation  is 


u.  .  -  c2u 
tt  xx 


.  P(t) 

m 


t  >  0, 

0  <  x  <  l 


(5) 


The  initial  conditions  and  boundary  conditions  for  the  free-free 
system  are 

I.C.'s:  u ( x  ,o )  =  o  ut(x,o)  =  o  (6) 

B.C.'s:  ux(o,t)  =  ^41  uv(£,t)  =  o  (7) 

The  solution  of  equations  (5)  through  (7)  gives  the  deformation  of 
any  point  for  any  time  greater  than  zero.  The  stress  is  then  simply 
obtained  from  Hooke's  Law,  a  =  Eux> 

The  solution  method  proceeds  by  transformation  of  the  given  system  to 
a  system  with  homogeneous  boundary  conditions.  This  is  achieved  by 
considering  a  solution  in  the  form 

u(x,t)  =  v (x ,  t }  +  U(x,t) .  (8) 

Boundary  condition  equations  (7)  become 

vx(0’t}  =  eF  '  Ux(0>t) 

(9) 

vxU,t)  -  -uxu,t) 

These  equations  become  homogeneous  when  U(x,t)  is  taken  as 

2 

U(x,t)  =  (x  -  ^  (10) 
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With  equations  (8)  and  (10),  governing  equations  (5)  through  (7) 
become 

vtt  •  c2yxx  '  f(x't)  01  > 


To  this  point  the  analysis  is  identical  to  that  presented  in  reference 

K  In  tha‘  ana,J,5's.  P(t)  was  approximated  as  a  half  sine  wave  Posin  n/a; 
here,  the  breech  force  is  approximated  by  the  Fourier  expansion,  ” 


P(t) 


I  1 

= a 

0  L  i 

k=l  k 


Ak  cos  +  B|<  sin 


(20) 


For  the  linear  system,  the  response  due  to  equation  (20)  is  obtained  by 
superposition  of  individual  terms. 

(i)  for  the  uniform  force  A  : 

o 

From  equation  (14),  f(x.t)  -  -  fa  .  Then  equation  (18)  gives 

fn  "  0  for  n-  The  initial  conditions  are  vn  =  o,  and 
vn  -  o  from  equation  (17),  and  hence  vn(t)  =  o.  Then  the 


stress  is 


o  =  Eux  =  EUX  =  0  -  £)  J° 


(21) 


This  is  the  stress  due  to  the  "rigid  body"  acceleration  for 
constant  force. 

(ii)  for  the  cosine  terms,  Ak  cos 


2kirt 

a 


From  equation  (14)  we  obtain, 

2 

f(x.t)  =  j(idkA  (it)2 

l  EA  \  a  J 


Then  equations  (17)  and  (18)  give 
d 


2c 1 
EA£ 


•  A.  cos 

K  a 


dt  vn(o)  =  0 


f  (t) 
nv  ' 


8k2l 

a2n2EA 


A.  cos  2-^ 

K  a 


(22) 

(23) 
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With  equations  (22)  and  (23),  the  solution  of  equation  (16) 

“n‘  } 


Ka\ 

y  (tl  =  — n _ £ 

V1'  2,  2 


2kTTt 

cos - cos 

a 


where 


and 


K  =  -  „ 

n  2  2 


2£ 


n  a  EA 


0) 

Yn  =  JL. 

2kTT 


Then  the  stress  associated  with  A.  cos  terms 


is 


^  =  Eu  =  E(v  +  U  ) 


V*  -2- 

C_i  unA 

n=i 


x  'X' 
A, 


.  _  2k  TT  t 

(yn  ~  1 )  [cos— -cos“ 


*  77 


(24) 


(25) 


(26) 


sin 


+  ( 1  -*>£».  a* 


rvrx 

t 

(27 ) 


The  latter  term  in  equation  (27)  is  the  contribution  associ¬ 
ated  with  the  "rigid  body"  acceleration, 
iii)  for  sine  terms,  B,  sin  • 

K  a 


Equation  (14)  gives 

x2 

f(x.t)  -  !  x  ■  it  1 2ki2 ' 

l~ur-  \— 

From  equations  (17)  and  (18), 


2c 

EA-f 


) 


A  me  2kTt 

A.  COS  - 

K  « 


3t  »„<»>  -  - 


4k£B, 


EAn^ira 


_  8k2l 


f  m  =  OK  1  n  2k7Tt 

'nlw  o'  4  B.  sin  - 

n  cxVeA  k 


(28) 

(29) 


,.N 


The  solution  of  equation  (16)  with  equations  (28)  and 
(29)  is 

2K_kB,a2  f 

T 


vn ( t )  - 


y2  -  2 
Tn 


Yn  -  1 


s,n  V  +  TTT 


,  ,  -  sin  f.-'o) 

»  (y‘  -  1)  “  '  ' 


where 

un 

Rn  =  (31) 

The  stress  associated  with  Bk  sin  is 

oo 


Again  the  latter  term  in  equation  (32)  results  from  the  "rigid 
body"  acceleration. 


The  stress  at  any  point  is  obtained  as  the  superposition  of 
equations  (21),  (27),  and  (32). 
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RESULTS  AND  DISCUSS ION 


Equations  (21),  (27),  and  (3?,  for  the  stress  „a<e  ^ 

on  the  NPS  -  IBM  360/67  digital  ccmnuter  for  specific  values  of  Ak  and 

Bk  obtained  by  a  Fourier  decomposition  of  the  given  breech  force.  This 

decomposition  is  done  by  the  IBM  subroutine  FORIT  on  the  computer.  The 

present  analysis  resolved  the  breech  force  into  26  sine  and  cosine  har- 
monies,  i.e., 

26 

P(t)  "  Ao  +  Z  (  Ak  cos  +  \  sin  |  (33) 


The  magnitude  of  the  harmonic  coefficients  Afc  and  Bk  depend  upon  the 

nature  of  the  breech  force.  For  this  analysis,  the  force  given  in  Fin. 

1  of  Ref.  2  was  used.  Approximation  of  that  force  by  equation  (33)  leads 
to 

1 - — — - - - - 

L_ 

\ 

Bk 

k 

*k 

bk 

0 

14.1680 

o 

14 

'.3164 

-.0186 

1 

-8.6082 

9. 6^38 

15 

-.1682 

-.0940 

2 

-3.2336 

-4.3103 

16 

-.0552 

.0146 

j 

1 .3139 

-.8107 

17 

-.1202 

.0983 

4 

-.1162 

.1077 

18 

-.2053 

.0899 

5 

.1089 

-.1257 

19 

-.1821 

.0827 

6 

- . 2036 

.3394 

1 

20 

-.2162 

.0815 

7 

-.3478 

-.0886 

21 

-.2723 

.0286 

8 

-.3001 

-.1664 

22 

-.2164 

-.0294 

9 

- .  2229 

-.3071 

23 

-.1599 

.0250 

10 

.1728 

-.2293 

24 

-.2146 

.0638 

11 

.2260 

.1067 

25 

-.2589 

-.0065 

12 

-.0425 

.3243 

26 

-.2162 

.0034 

13 

-.3157 

.2287 
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A  plot  of  equation  (33)  with  the  above  values  of  and  Bk  is 
given  in  Fig.  1  for  a  .026  second  time  dura1"1''*  of  breech  force. 

Figures  2  through  5  are  the  results  for  the  stress  in  the  bar  due  to 
the  breech  force  shown  in  Fig.  1  as  a  function  of  time  for  j  ratios 
of  0,  1/4,  1/2,  and  3/4,  r°spectively.  (The  analysis  is  for  the  5  in. 
steel  shell  of  0.19  in.  thickness).  Comparing  Fig.  l  with  Figs.  2 
through  5  reveals  that  the  internal  force  appears  as  a  duplicate  of 
the  breech  force,  with  the  stress  curves  scaled  versions  of  the  force 
input.  In  Figs.  6  through  8,  the  contributing  harmonics  are  individu¬ 
ally  displ  ed.  The  elastic  dynamic  effects  are  the  small  oscillations 
about  each  harmonic  itself. 

Figures  !0  through  14  are  the  results  of  the  analysis  for  a 
*ime  duration  of  0.0052  seconds  for  the  brech  force.  This  is  a  re¬ 
duction  in  the  actual  duration  by  a  factor  of  five.  This  sample  prob¬ 
lem  was  considered  in  order  to  assess  the  influence  of  the  magnitude 
of  the  force  duration.  The  elastic  dynamic  effects  can  be  clearly 
seen  in  these  figures.  Note  that  the  magnitude  of  stress  has  beer, 
significantly  changed  by  the  dynamic  effects. 

CONCLUSIONS 

The  results  of  this  analysis  (Figs.  1-5)  show  that  dynamic 
effects  are  negligible  for  the  breech  force  of  0.026  seconds  time 
duration  of  Figure  1,  Reference  2.  Hence,  when  the  propellant  is 
not  bonded  to  the  motor  case,  the  case  responds  to  the  breech  force 
as  if  it  were  applied  in  a  static  sense,  i.e.,  the  maximum  internal 
axial  load  in  the  case  can  be  computed  on  the  basis  of  the  mass 
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distribution  and  the  maximum  rigid  body  acceleration  of  the  rocket. 
This  is  the  same  conclusion  arrived  at  in  Ref.  1  based  on  a  simple 
one-half  sine  wave  pulse;  hence,  accounting  for  the  more  complicated 
actual  breech  pressure  pulse  confirmed  the  original  conclusion. 

Further  investigation  (Figs.  10-14)  verified  the  dependence 
of  the  dynamic  effects  on  the  presence  of  significant  higher  harmonic, 
terms  in  the  breech  force.  Such  harmonics  may  be  the  result  of  shor¬ 
tening  the  time  duration  cf  loading,  or  by  the  presence  of  a  local 
rapid  perturbation  of  force  on  an  otherwise  acceptable  smooth  breech 


force  input. 
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Figure  5 
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Figure  10 
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Figure  11 
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Figure  13 
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Figure  14 
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